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Abstract: We present a complete method to construct QCD-protected observables based 
on the exclusive 4-body B-meson decay B^ — )• K*^i~^£~ in the low dilepton mass region. 
The core of the method is the requirement that the constructed quantities should fulfil the 
symmetries of the angular distribution. We have identified all symmetries of the angular 
distribution in the limit of massless leptons and explore: a new non-trivial relation between 
the coefficients of the angular distribution, the possibility to fully solve the system for the 
K* amplitudes and the construction of non-trivial observables. 

We also present a phenomenological analysis of the new physics sensitivity of angular 
observables in the decay based on QCD factorisation. We further analyse the CP-conserving 
observables, A!^ , A^ and A!^ . They are practically free of theoretical uncertainties due 
to the soft form factors for the full range of dilepton masses rather than just at a single 
point as for ^fb- They also have a higher sensitivity to specific new physics scenarios 
compared to observables such as ^fb- Moreover, we critically examine the new physics 
reach of CP-violating observables via a complete error analysis due to scale dependences, 
form factors and vlQco/^T-fo corrections. We have developed an ensemble method to evaluate 
the error on observables from AqcB/'nT'b corrections. Finally, we explore the experimental 
prospects of CP-violating observables and ffiid that they are rather limited. Indeed, the 
CP-conserving (averaged) observables A^ (with i = 2,3,4) will offer a better sensitivity 
to large CP phases and may be more suitable for experimental analysis. 

Keywords: B-Physics, Rare Decays. 



Contents 

1. Introduction 2 

2. Theoretical framework 4 

2.1 Differential decay distribution 4 

2.2 QCDf/SCET framework 6 

2.3 Estimating ^Qco/'^T-fe corrections 8 

3. Symmetries and observables 9 

3.1 Infinitesimal symmetries 10 

3.2 Explicit form of symmetries 11 

3.3 Relationship between coefficients in differential distribution 11 

3.4 Experimental issues 13 

3.5 Constructing observables 13 

3.6 More general cases 14 

4. Experimental Sensitivities 15 

4.1 Experimental Analysis 16 

4.1.1 Generation 17 

4.1.2 Observable sensitivities 17 

4.1.3 CP asymmetries 17 

4.2 The polynomial ansatz re-examined 18 

4.3 Fit Quality 18 

4.4 Discussion 20 

5. Analysis of CP-violating observables 20 

5.1 Preliminaries 21 

5.2 Phenomenological analysis 22 

6. Analysis of CP-conserving observables 27 

6.1 Leading-order expressions of ^^ 27 

(5) 

6.2 Leading-order expressions of ^^ 29 

6.3 Analysis of A^ and a':^^ 31 

7. Conclusion 31 

A. Kinematics 33 

B. Theoretical input parameters and uncertainties 34 



1. Introduction 

The LHC era is just beginning. Flavour physics will play an important complementary role 
to direct searches for the theory that lies beyond the standard model (SM). One central 
strategy in this period is to construct observables that are mostly sensitive to specific 
types of new physics (NP), in such a way that a deviation could immediately provide 
information on the type of NP required: isospin breaking NP, presence of right-handed 
currents, scalars, etc. It is essential to work in a bottom up approach in the direction of 
constructing a decision tree that help us to discern which features the NP model must 
incorporate and then try to match them into a group of models. 

Few decays are able to provide such a wealth of information with different observables 
as Bd — )• K*^t^i.~ , ranging from forward-backward asymmetries (j4fb) and isospin asym- 
metries to a large number of angular observables. Each of these observables can provide 
information on the different types of NP mentioned above. First published results from 
BELLE [1] and BABAR [2] based on O(IOO) decays already demonstrate their feasibility. 

In the early years of LHC running one will be restricted to those observables that 
may be extracted from the angular distribution using relatively simple analyses. A study 
of those observables relevant for the first few fb~ may be found in [3]. However, once 
enough statistics have been accumulated to perform a full angular analysis based on the full 
4-body decay distribution of the B^ — >• K*^(.'^(.~ , one has the freedom to design observables 
with reduced theoretical uncertainties and specific NP sensitivity. 

In [4], it was proposed to construct observables that maximise the sensitivity to con- 
tributions driven by the electro-magnetic dipole operator Oj, while, at the same time, 
minimising the dependence on the poorly known soft form factors. This led to the con- 
struction of the observable A)^ , based on the parallel and perpendicular spin amplitudes 
of the K*^ . The basic idea behind the construction of the observable was inspired by the 
zero point of A-pB when calculated as a function of the dilepton mass squared, q^ . The 
zero point has attracted a lot of attention because of its cleanliness; only at that point one 
gets a complete cancellation at LO of the form factor dependence and its precise position 
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may provide information on the fundamental theory that lies beyond the SM. For A)j, the 
soft form factor dependence cancels at LO, not only at one point, but in the full cf' region 
thus providing much more experimental information. Moreover, the angular observable is 
highly sensitive to new right-handed currents driven by the operator Oy [5], to which ^fb 
is blind. 

Looking for the complete set of angular observables sensitive to right-handed currents, 
one is guided to the construction of the so-called As^ and A^ which include longitudinal 
spin amplitudes [6]. The observables Aj^ (with i = 2,3,4) use the K*^ spin amplitudes as 
the fundamental building block. This provides more freedom to disentangle the information 
on specific Wilson coefficients than just restricting oneself to use the coefficients of the 
angular distribution as it was recently done in [7]. For instance, Ar^\ being directly 
proportional to Cy enhances its sensitivity to the type of NP entering this coefficient. 
Moreover, using each coefficient of the angular distribution instead of selected ratios of 
them induces a larger sensitivity to the soft form factors. 



The spin amplitudes are not directly observable quantities; to ensure that a quantity 
constructed out of the spin amplitudes can be observed, it is necessary that it fulfils the 
same symmetries as the angular distribution. This observation has the important conse- 
quence [6] that A>p (first proposed in [8]) cannot be extracted from the angular distribution 
because it does not respect all its symmetries. Only a measurement of definite helicity dis- 
tributions would allow it, but that is beyond any particle physics experiment that can 
currently be imagined [6]. 

To identify all the symmetries of the angular distribution is one of the main results of 
this paper. We discuss the counting of all the symmetries of the distribution in different 
scenarios, with and without scalars and with and without mass terms. We explain the 
general method of infinitesimal transformations that allow us to identify all the symmetries, 
and we develop here in full detail the explicit form of the four symmetries in the massless 
case with no scalars. As an important cross check of this result, we solve explicitly the set 
of spin amplitudes in terms of the coefficients of the distribution, making use of three out 
of the four symmetries. Two important consequences of this analysis are: in solving the 
system one naturally encounters an extra freedom to fix one of the variables, and there is a 
non-trivial constraint between the coefficients of the angular distribution considered before 
as free parameters. It is remarkable that this unexpected constraint is valid for any decay 
that has this same structure. 

Finally, we provide an illustrative example of the use of the method of designing 
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observables with an observable called A^ that mixes simultaneously left/right and per- 
pendicular/parallel spin amplitudes in a specific way that none of the coefficients of the 
angular distribution exhibits, opening different sensitivities to Wilson coefficients. 

In the second part of the paper we present a phenomenological analysis of the various 
angular observables based on a QCD factorisation (QCDf) calculation to NLO precision. 
Recently, a very detailed analysis of angular quantities of the decay Bfi — t- K*^h'^ij,~ 
in various NP scenarios [7] and also an analysis of the NP sensitivities of angular CP 
asymmetries [9] were presented. In contrast to the former work [7], we do not assume that 
the main part of the vIqcd/"^;, corrections are inside the QCD form factors, but use the soft 
form factors and develop a new ensemble method for treating these unknown corrections in 
a systematic way. The main differences to the latter analysis of CP violating observables is 
the redefinition of the CP asymmetries in order to eliminate the soft form factor dependence 
at LO and the inclusion of the Aqcd/tjii, corrections into the error budget, which turn out 
to be significant in the presence of new weak phases. 

In [6] the experimental preparations for an indirect NP search using these angular 
observables were worked out, showing that a full angular analysis of the decay B^ — )• 
K*^fi'^fi~ at the LIIC6 experiment offers great opportunities. We re-evaluate this analysis 
in light of the fourth symmetry for the angular distribution and conclude that it has no 
effect on the estimated experimental errors as all observables are indeed invariant under 
this symmetry. We extend the experimental sensitivity study to CP-violating observables 
and show that even with an upgraded LHC6 there is no real sensitivity to CP-violating 
NP phases in Cg and Ciq. 

The paper is organised as follows: Sec. 2 briefly recall the differential distribution 



in Bd — >■ K*^i'^i^ and the theoretical framework of QCDf and soft-colhnear effective 
theory (SCET), Sec. 3 extends and completes our previous discussion about symmetries 
in the angular distribution, its experimental consequences are discussed in Sec. 4, and we 
perform a phenomenological analysis of the CP-violating and CP-conserving observables 
in Sees. 5 and 6 respectively. 

2. Theoretical framework 

The separation of NP effects and hadronic uncertainties is the key issue when using flavour 
observables in a NP search. Our analysis is based on QCDf and SCET and critically 
examines the NP reach of those observables via a detailed error analysis including the 
impact of the unknown ^QCD/"^b corrections. In order to make the paper self contained, 
we briefly recall the various theoretical ingredients of our analysis. 

2.1 Differential decay distribution 

The decay B^ — t- K*^i~^i~ , with K*^ — t- K^-k^ on the mass shell, is completely described 
by four independent kinematic variables, the lepton-pair invariant mass squared, g^, and the 
three angles 9i, 6k, 4>- Summing over the spins of the final state particles, the differential 
decay distribution of B^ — >■ K*^i~^£~ can be written as 

:r2-, — f r fl ^A = ^-^^i'^ ^'' ^^' '^) ' (2.1) 

ag^ a cos ffi a cos Ok dcp Szvr 
The dependence on the three angles can be made more explicit: 

j{q\eueK,(t>) = 

= Jis siv? Ok + Jic cos^ Ok + ( ^2* sin^ Ok + J2c cos^ Ok) cos 20i + J3 sin^ Ok sin^ Oi cos 2(/) 
+ J4 sin 20 k sin 20i cos (/) + J5 sin 20k sin 9i cos (p + {Jqs sin^ Ok + Jgc cos^ Ok) cos Oi 

+ J7 sin 20k sin Oi sin </> + Jg sin 20k sin 20i sin </> + Jg sin^ Ok sin^ 0i sin 2<p ■ (2.2) 

As the signs of the expression depend on the exact definition of the angles, we have made 
their definition explicit in Appendix A. 

The Ji depend on products of the six complex K* spin amplitudes, Au ' , Aj^ and 

Aq' in the case of the SM with massless leptons. Each of these is a function of q . The 
amplitudes are just linear combinations of the well-known helicity amplitudes describing 
the B —7- Ktt transition: 



^^,11 = {H+i T H^i)/V2 , Ao = Ho. (2.3) 

Two generalisations will be made from the massless case within our analysis: if the leptons 
are considered massive the additional amplitude At has to be introduced. And if we allow 
for scalar operators, there is a new amplitude As- Both can be introduced independently 
of the other. For the Jj we find the following expressions (see also [10, 11, 12, 4]) *: 

,^^a= ^l±l!l \\Ai\^ + |A[|2 + (L ^ R)] + ^Re f^i^f + A^^Af) , (2.4a) 
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'The generalizations to the case which includes scalar operators was recently presented in [7] 
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(2.4c) 
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(2.41) 

(2.5) 



The notations with the letters a-m has been included to make the comparison to [6] easier. 
Note that Jqc = in the massless case. 

The amplitudes themselves can be parametrised in terms of the seven B — t- K* form 
factors by means of a narrow- width approximation. They also depend on the short-distance 
Wilson coefficients Ci corresponding to the various operators of the effective electroweak 
Hamiltonian. The precise definitions of the form factors and of the effective operators are 
given in [6]. Assuming only the three most important SM operators for this decay mode, 
namely Oy, Og, and Oio, and the chirally fiipped ones, being numerically relevant, we 
have ^ 
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^Following common convention, we use the effective Wilson coefficients of these operators which include 
contributions from four-quark operators as well. 
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where the Ci denote the corresponding Wilson coefficients and 
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Finally we note that, if one additionally considers scalar operators then At is modified 
by the new Wilson coefficients and an additional amplitude, As, proportional to the form 
factor Ao('?^)) is introduced. 

2.2 QCDf/SCET framework 

The up-to-date predictions of exclusive modes are based on QCDf and its quantum field 
theoretical formulation, soft-collinear effective theory (SCET) [13, 14]. The crucial the- 
oretical observation is that in the limit where the initial hadron is heavy and the final 
meson has a large energy [15] the hadronic form factors can be expanded in the small 
ratios ^Qco/^ife and Aqcb/E, where E is the energy of the meson that picks up the s 
quark from the B^ decay. Neglecting corrections of order 1/mi, and Og, the seven a-priori 
independent B — )• K* form factors reduce to two universal form factors .^_l and ^u [15, 16]. 
These relations can be strictly derived within the QCDf and SCET approach and lead to 
simple factorisation formula for the B — t- K* form factors 



Fi(g2) =H^^ + <^B®Ti0<^K*+ 0{AqcB/mb) 



(2.9) 



There is also a similar factorisation formula for the decay amplitudes. The rationale of such 
formulae is that the hard vertex renormalisations (Hi) and the hard scattering kernels (Tj) 
are quantities that can be computed perturbatively so they can be separated from the non- 
perturbative functions that go with them; i.e. the light-cone wave functions ($i) which are 
process-independent and the soft form factors (^) which enter in several different B — ?• K* 
processes. 

In general we have no means to calculate AqcD/iTT-b corrections to the QCDf amplitudes 
so they are treated as unknown corrections, with the method used for this described in the 
following section. This, in general, leads to a large uncertainty of theoretical predictions 



based on the QCDf/SCET which we will explore systematically and make manifest in our 
phenomenological analysis. 

We do not follow here the approach of [7] where the full QCD form factors are used in 
the QCDf formulae. There it is assumed that the main part of the ^QCD/'"^fe corrections 
are inside the QCD form factors, and additional ^QCD/"^fe corrections are just neglected. 
Clearly some of the AQQY)/miy corrections could be moved into the full QCD form factors. 
However, there is no robust quantitative estimate of the additional corrections and, thus, 
it is not allowed to neglect those unknown corrections, especially in view of the expected 
smallness of new physics effects. 

We follow here another strategy. We construct observables in which the soft form 
factor dependence cancels out at leading order. Then the influence of the soft form factors 
to the physics is almost eliminated from the phenomenological analysis in a controlled way. 
On the other hand we make the uncertainty due to AQCB/iTT-b corrections manifest in our 
analysis. It is not expected that there are as large as 20 — 30% as in the B — )• tttt decay 
as argued below. The inclusion of the 5 — 10% errors due to the ^qcd/?^6 corrections 
in our analysis is exploratory of its impact on our observables, even at the risk to be too 
conservative. Obviously, it is this issue which calls for improvement in view of the new 
physics reach of these modes. 

The theoretical simplifications of the QCDf/SCET approach are restricted to the kine- 
matic region in which the energy of the K* is of the order of the heavy quark mass, 
i.e. q'^ <^ m^. Moreover, the influences of very light resonances below 1 GeV^ question the 
QCDf results in that region. In addition, the longitudinal amplitude in the QCDf/SCET 
approach generates a logarithmic divergence in the limit (7^ — t- indicating problems in 
the theoretical description below 1 GeV^ [13]. Thus, we will confine our analysis of all 
observables to the dilepton mass in the range, 1 GeV^ ^ g^ ^ 6 GeV^. 

Using the discussed simplifications, the K* spin amplitudes at leading order in l/nif, 
and Os have the very simple forms: 
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(2.10b) 

C\\{Ek*), 
(2.10c) 

(2.10d) 
The scalar spin 



with s = q /m^, rhi = rrii/mB- Here we neglected terms of 0{mj^t] 
amplitude As is also proportional to ^\\{Ek*) in this limit. 

The symmetry breaking corrections of order q^ can be calculated in the QCDf/SCET 
approach. Those NLO corrections are included in our numerical analysis following [13, 14]. 



They are presented in the Appendix of [6]. 



2.3 Estimating ^qcd/'ti6 corrections 

Our observables have reduced theoretical uncertainties due to the cancehation of the soft 
form factors. However, the relations used to make these cancellations are only valid at 
LO in the ^qcd/^?) expansion, and corrections to higher orders are unknown. For these 
theoretically clean observables to be useful, the impact of these corrections on the observ- 
ables must be robustly bounded. If NP is to be discovered in Bd — )• K*^i^i~, it must 
be possible to demonstrate that any effect seen is indeed NP and not just the effect of an 
unknown SM correction. 

To evaluate the effect of the ^qcd/'"^6 corrections, we parametrise each of the K*^ 
spin-amplitudes with some unknown linear correction. 



A', = M1 + Cie 



idi 



(2.11) 



where Ci is the relative amplitude and 6i the relative strong phase. If we vary Cj and 
9i within their allowed ranges, an estimate for the theoretical uncertainty due to these 
unknown parameters can be found. In order to make this parametrisation generic, however, 
extra terms must be introduced. In principle the effective Hamiltonian which controls the 
decay has three terms. 



^cff = n 



(u)SM 
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(t)SM 



+ n 



(t)NP 
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The first term is very small as it is suppressed by the factor A^ = KibK*s/^tbKs but is 
responsible for all the SM CP-violation in the decay; the second term is responsible for 
the decay in the SM; and the third adds possible NP contributions. A fourth possible 



term % 



(«)NP 
off 



generically does not contribute to the model independent amplitudes and is 



neglected. Each of these contributions is generated by different sets of diagrams and may 
have different values of Ci and 6i. 

Each amplitude must be modified to include the three sub-amplitudes with their cor- 
rections: 
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(2.13) 



It is assumed that only a single NP operator is active so as not to introduce extra terms. In 
this formalism, the SM CP-violating, SM CP-conserving, and NP parts of the amplitude 
are then allowed to have independent AQQD/rrH) corrections and strong phases. 

An estimate of the theoretical uncertainty arising from the unknown ^qcd/"^6 correc- 
tions and strong phases can now be made using a randomly selected ensemble. For each 
member of the ensemble, values of Ci_3 and ^1-3 are chosen in the ranges Ci G [—0.1, 0.1] 
or Ci G [—0.05,0.05] and 6i G [— 7r,7r] from a random uniform distribution. This is done 
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for the seven amplitudes. At, Aq' , A..' , ^ ,' , to provide a complete description of the 



decay. It is assumed that the corrections and phases are not functions of q^, although in 
practise they may actually be. Any unknown correlations are also ignored. While these 
effects could lead to an underestimate of the theoretical envelope, it is thought that this 
method allows for a conservative estimate of the theoretical uncertainties to be made. 

To estimate the contribution to the theoretical uncertainties from ^qcd/'"^6 corrections 
for a particular observable, each element in the ensemble was used to calculate the value 
of that observable at a fixed value of g^. A one a error is evaluated as the interval that 
contains 66% of the values around the median. This is done for both d E [—0.05, 0.05] and 
Ci G [—0.1,0.1] to illustrate the effects of five and ten percent corrections. By repeating 
this process for different values of (f' , bands can be built up. No assumption of Gaussian 
statistics has been made; the bands illustrate the probable range for the true value of each 
observable, given the current central value. The method allows for the probability that a 
given experimental result is due to an unknown SM correction to be found. 

The choice |Cj| < 10% is based on a simple dimensional estimate. We emphasize here 
that there is no strict argument available to bound the AQQ^/rrn, corrections this way. But 
we can state that the chiral enhancement of ^qcD/'^b corrections in the case of hadronic 
B decays does not happen in the case of the semileptonic decay mode with a vector final 
state. 

The process described here avoids any assumptions about correlations between the 
corrections and is thus statistically more rigorous than what was done in [6], where cor- 
rections to amplitudes were considered one by one and then added in quadrature. The 
^QCD/'iT'b bands it produces are reduced when compared to those of [6]. It also allows us 
to investigate the effect of the AQQ^/mi^ corrections for CP-violating observables. 

3. Symmetries and observables 

The experimental degrees of freedom determined by the J, terms and the theoretical degrees 
of freedom determined by the spin amplitudes Aj have to match. There are two effects to 
consider for this: different values of the Aj can give rise to the same differential distribution 
Eq. (2.1) and thus cannot be distinguished; and in some cases the experimental coefficients 
are not independent, meaning that not all arbitrary values of the Ji are possible. The first 
effect we call a continuous symmetry transformation. For the degrees of freedom to match 
we have 

nc-nd = 2nA - rig , (3.1) 

where ric is the number of coefficients in the differential distribution (the number of Jj), 
rid the number of dependencies between the different coefficients, ua the number of spin 
amplitudes (the Aj, each is complex and hence has two degrees of freedom), and rig the 
number of continuous symmetries. 

We considered this situation in our previous paper [6] for the case of massless leptons 
and return to it again here. It is easy to see that in the massless limit, Jis = 3J2s and 
Jic = —J2c- What is not so obvious is that Jg can be expressed in terms of the other 8 
remaining coefficients. Going back to Eq. (3.1) it can be seen that the massless case in fact 
must have 4 symmetries and not 3 as we claimed in the previous paper. 



Below we first outline how the symmetries and dependencies can be identified before 
we move onto their explicit form and the interpretation. 

3.1 Infinitesimal symmetries 

By an infinitesimal symmetry is meant one where the theoretical spin amplitudes Aj are 
changed in an infinitesimal way leaving the Jj coefficients in Eq. (2.4) unchanged. The 
infinitesimal symmetries will define a system of coupled ordinary differential equations 
that, if solved, are the global symmetries we look for. There is no guarantee that these 
symmetries will allow for the continuous transformation between two arbitrary sets of 
amplitudes which have the identical angular distribution; there could in principle be several 
disjoint regions separated by divergences. 

If we, in this example, look at massless leptons and ignore the scalar amplitude, the 
coefficients of the spin amplitudes can be defined as a vector A with 12 components 



A = (Re{A^),lm{A{),Re{A^),lm{A^),Re{A^),lm{A^), 

Re(A^), Im(^^), Re(^f ), Im(^f ), Re(^^), Im(A^)) (3.2) 



corresponding to the real and imaginary parts of the amplitudes. For each of the coefficients 
Jj we can find the derivative with respect to the spin amplitudes. As an example 

V( Jic) = (0, 0, 0, 0, 2Re(^^), 2Im(yl^), 0, 0, 0, 0, 2Re(<), 2Im(<)) . (3.3) 

There will be eleven such gradient vectors in the massless case, as Jqc = 0. 
Now, any infinitesimal transformation can be written on the form 

A' = A + 6s. (3.4) 

For the infinitesimal transformation to leave the coefficients unchanged, the vector 5s has 
to be perpendicular to the hyperplane spanned by the set of gradient vectors. Or in other 
words, 5s represents a symmetry if, and only if 

WiG Ji:Vi± 5s. (3.5) 

Looking back at Eq. (3.1) we have, for the massless case, Uc = 11. If the Jj were 
all independent the gradient vectors would span an 11 dimensional hyperplane. In fact, it 
turns out that they only span 8 dimensions"'" , which shows that there are three dependencies 
between the Jj's, giving n^ = 3. As we have ha = 6 from the amplitudes we see from 
Eq. (3.1) that we have n<( = 4 corresponding to 4 symmetries. For the dependencies, only 
the first two Jis = 3J2s and Jic = —J2c are trivial; the third one we derive in the next 
section. 



"'"Any program able to handle symbolic algebra will be able to show this. 
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3.2 Explicit form of symmetries 

It is helpful for the discussion to make the following definitions. 

n2 = (^i,-^f), 



or in terms of helicity amplitudes 



mi 



771,2 



1 
1 

7^ 



R *•] 



{ni + n2) = {H^i,H^i ), 



{ni-n2) = {H^i,H^i ), 



ms = n3 = iH^,H^*) 



(3.6a) 
(3.6b) 
(3.6c) 



(3.7a) 

(3.7b) 

(3.7c) 



In fact, all the information of the angular distribution is encoded in the moduli of the three 
77j vectors and their relative complex scalar products: 



ni • 772 



3 

Jqs 



Jls — Js , 
iJg, 



P2I 

771 -ris 



- Jls + J3 , 



V2J4-i 



J7_ 

V2 



I^-Sl =Jlc; 
772 • 773 = 



(3.8) 
^-iV2Js, (3.9) 



where 77^ being a complex vector implies that the scalar product is 77j • 77j = J2k^ik''^*j ■ 
The coefficients J2s and J2c are absent because they are obviously redundant. 

The differential distribution is invariant under the following four independent symme- 
try transformations of the amplitudes 



n,- 



e"PL 
e" 



cos( 
sin( 



- sm( 
cos( 



cosh iO — sinh i6 
- sinh id cosh ii 



Ui 



(3.10) 



where ^i,^ </>/?, and Q can be varied independently. Identical transformations can be 
carried out on the 777^. Normally, there is the freedom to pick a single global phase, but as 
L and R amplitudes do not interfere here, two phases can be chosen arbitrarily as reflected 
in the first transformation matrix. 

The interpretation of the third and fourth symmetry is that they transform a helicity 
+1 final state with a left handed current into a helicity —1 state with a right handed 
current. As we experimentally cannot measure the simultaneous change of helicity and 
handedness of the current, these transformations turn into symmetries for the differential 
decay rate. 

3.3 Relationship between coefficients in differential distribution 

As was mentioned earlier, we have identified an extra dependency among the coefficients 
in the massless case. Here we outline how it can be derived. 
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If we use the two global phase symmetry transformations we can rotate the vector ni 
to make it real {Au and Af^ become real).§ We can then choose the angle 6 of the third 
symmetry to make Au = 0. Notice that we have not made use of the fourth symmetry. 
The implications of this fourth symmetry will become manifest when solving the system. 
With these choices 

ni = {0,Af), (3.11) 

where A^ is a positive real parameter. Using three of Eqs. (3.8)-(3.9) together with the 
symmetries, one can determine four of the spin amplitudes (their moduli and phases): 

AI( = 0, (3.12a) 



Af = ^ArH\' = \|lJls-J3, (3.12b) 

A5 = -^^ = -i^^?^^Sl, (3.12c) 

A^ = ^,= '''-''' (3.12d) 

The remaining three equations from Eqs. (3.8)-(3.9) determine, on one side, the moduli of 
A^ and A^, 



ii|2_i„_i2 |( ni.re2)| ^ 9 Jig - Ji, - iJqs ~ Jg 
Fl| 3'^^s — Js 



Ai\^ = |n2p - '^'7"7/ = 9-1- ;3 4-6. -9 ^ (3^^3^) 



K ni.ns)^ ^ Jic (i Jis - Js) - 2J | - ij| 

Pi I S-^ls ~ "^3 



1^0 I — PsI nr~|2 — ~ 2~;^ 7 ' (3.13bj 



and on the other, the phase difference corresponding to the previous two amplitudes: 



f,L) _ (n2 • ra3)|?7-i|^ - (ra2 • ni)(ni • ^3) 



[(|ni|2|n2P - |(n2 • ni)|2) (|ni|2|n3|2 - |(n3 • ni)p)]^''^ 
_ "^5 (3'^is ~ "^3) — JaJqs — JjJg — i {-^JisJs — 2J3 Js + 2J4J9 — 2^63^7) 

[2 (g-^L - Ji - i-^ls - J9) {Jic {^Jis - J3) - 2J| - i^Jr)] 

(3.14) 

Here is where the fourth symmetry becomes manifest. On one side, this equation tells 
us that you have the freedom to choose one of the two phases cp^ or (pQ to zero. On the 
other side, given that the LHS of the previous equation is a pure phase, the modulus of the 
RHS should be one. This implies the following important non-trivial relationship between 
the coefficients of the distribution 

(2 Ji, + 3 J3) (4 J| + J|) + (2 Ji, - 3 J3) {Jj + 4 J|) 
'' 16 Jf 2 - 9 (4 J32 + Ji^ + 4 J|) 

16 J,2, - 9 (4 J| + J|, + 4 J|) • ^'-^'^ 



^Indeed the system can also be solved using only one of the two global symmetries and keep A^ complex. 
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It is important to remark that this equation can be very easily generahsed to the massless 
case with scalars using the relations Ju = 3J2s and Jic = —J2c in the previous equation. 
Also the massive case with no scalars can be included by introducing the /3 factors inside 
the Jj coefficients. There is no such equation in the massive case with scalars due the fact 
that the number of coefficients of the experimental distribution is identical to the number 
of theoretical amplitudes and symmetries (see Tab. 1). 

3.4 Experimental issues 

The symmetries discussed above can be used to fix the spin-amplitude components by 
choosing specific values of the relevant rotation angles. We give an explicit example of 
this for the case where the lepton mass is neglected. We choose to make the following 
constraint: 

Re(A{') = Im(^f ) = lm{Af) = Im(^i) = 0. (3.16) 

This can be achieved by first performing the last transformation, shown in Eq. (3.10), 
with the value of 9 given by: 



sm0=A/^, ^°^^~=V^' ^^-^^^ 



where 



A 



1 + 4 



Re(A|^)Im(A|f ) + Re(A^)Im(A|^) 



Re{Afy + Im{Afy - Re(^f )2 - Im(ylf )^ 
Next, the third rotation angle, 0, is used again in Eq. (3.10): 



(3.18) 



yrT^Re(^[) - \/z^Ilm(A^) 
tan9=^^= — 1^- ; -4t. (3.19) 

The L- fields are phase-shifted by (p^: 

cos (9[cos e lin{A\ ) - sin 6* Im(^f )] + sin 6i[cos 9 ReiA^) + sin 9 Re(^f )] 

tan (pL = = 1 r ~ r ir~ ; 

sm.9[cos9lm{Al) - sin6' Im(^^)] + cos6'[cos6' Re(^^) + sm9Re{Af)\ 

(3.20) 
and finally the last i?-field transformation can be performed substituting (-L— s- ||) and 
(L o R) into the previous expression: 

cos^[cos6' Im(A?) - sin6' Im(yl|f)] + sin6i[cos6' Re(A|f ) + sin6' Re(^?)] 

tan cbR = ~ — . 

sin6l[cos6' Im(^{') - sin 6* Im(^p] + cos [cos 6* Re (yip + sin 6* Re(y4|f )] 

(3.21) 

3.5 Constructing observables 

In [4, 6], as well as here, we use the spin amplitudes to construct our observables. There 
are two main advantages of this approach, one is experimental and the other is theoretical. 
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On the experimental side, we have found that fitting directly the angular coefficients Jj, 
without taking into account the relations between them, leads to fit instabilities. These 
relations, coming from the underlying K*^ spin amplitudes, can be found in Sec. 3.3. 
The theoretical argument has to do with our aim at constructing observables that fulfil 
certain criteria, namely maximal sensitivity to a specific NP operator, like new right-handed 
currents, and minimal sensitivity to poorly known form factors. Given that our main tools 
are the spin amplitudes, it is a straight-forward exercise to design observables with a specific 
NP sensitivity and small hadronic uncertainties. We also have more freedom to construct 
observables than just using each coefficient of the distribution as an observable. As the 
spin amplitudes can be extracted directly in the full-angular analysis, there is no penalty 
on the final experimental uncertainty from using a non-trivial functional form to make the 
observable. 

The symmetries of the angular distribution play a crucial role in our approach. Once a 
quantity has been designed, it is a necessary condition for being an observable based on the 
angular distribution that it respects all the symmetries of this distribution. For example 
in [6] , we have shown explicitly that a previously discussed transversity amplitude A>p does 
not fulfil all the symmetries of the angular distribution. This implies that this quantity 
cannot be measured at the LHCb experiment or at future super- i? factory experiments; a 
measurement of the spins of the final-state particles would be required for that. 

Let us finally discuss a new CP-conserving observable that we call A)j. . It is defined 
as: 

(5^ \A^Af + Af*Af\ 

T — 77772 I ,p|2 7772 7772 • (3.22) 

l^il +l^±l +I^S'I +l^fl 

(2) 

It probes the transverse spin amplitudes A± and An in a different way than A}^ . Direct 
inspection of Eq. (2.4) shows that there is no single angular coefficient mixing L with R 

(5) 

and _L with || simultaneously in the way ^r^ does. 

It is a simple exercise to check that this observable fulfils the four symmetries described 
in Eq. (3.10). Once this invariance is fulfilled" one is allowed to use the explicit solution 
in the massless case provided in the previous subsection Eqs. (3.8)-(3.9): 



.(5) 



^ Vl6Jf -9J|^-36(J| + J^) _ ^^23) 

A discussion on the properties and sensitivities of this observable is presented in Sec. 6. 

3.6 More general cases 

The discussion of the differential symmetries from Sec. 3.1 can be generalised to the cases 
where the leptons are no longer considered massless and where a scalar amplitude is in- 
cluded: 



^Notice one could try to write Al^ in terms of the Ji using the expUcit solution, but this is not allowed 
since AJj, is not invariant [6] . 
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Massless leptons with scalars The inclusion of the scalar amplitude As, gives us seven 
amplitudes. The four explicit symmetries in Eq. (3.10) are still valid and we have in 
addition 

A's = e^'t''As, (3.24) 

expressing that the phase of As cannot be determined. 

Massive leptons without scalars We have the seven amplitudes A_^ , Au ' , Aq' and 
At in this case and still eleven coefficients. As a fact of elementary quantum mechanics 
we still have a global phase transformation corresponding to (j)^ = (j)R, but the other 
two symmetries from the massless case are no longer valid. There is a new symmetry 
concerning the phase of At given as: 

A[ = e'^'At. (3.25) 

This leaves us with two symmetries where only the differential form is known. 

Massive leptons with scalars We now have all eight amplitudes and, with the inclusion 
of Jqc, we have twelve coefficients. The global phase transformation, 0^ = </)/?, and 
the phase transformation of At in Eq. (3.25) are still valid. In this case, there is no 
dependency between any of the coefficients, leaving us with two symmetries where 
only the differential form is known. 

So while we in some cases only know the differential form of the symmetries, we are still 
able to test if observables respect the symmetries (see Sec 3.5) and we can also determine 
the optimal set of amplitudes to fit for in an experimental fit (see Sec. 3.4). In Tab. 1 we 
summarise the full knowledge about the symmetries. 



Case 




Coefficients 


Dependencies 


Amplitudes 


Symmetries 


rui = 0, As = 


= 


11 


3 


6 


4 


m£ = 




11 


2 


7 


5 


me >0, As = 


= 


11 


1 


7 


4 


■me> 




12 





8 


4 



Table 1: The dependencies between the coefficients in the differential distribution and the sym- 
metries between the amplitudes in several special cases. 



4. Experimental Sensitivities 

In [6], a fitting technique was investigated that allowed the extraction of the K*^ spin 
amplitudes from the full angular distribution in the massless lepton limit. Eq. (2.1) can 
be interpreted as a probability density function (PDF) and normalised numerically. We 
parametrise it in terms of six complex K*^ spin amplitudes, which are functions of q'^ only. 
In the limit of infinite experimental data, and for a fixed value of q'^, these amplitudes can 
be found by fitting the relative contribution of each angular coefficient as a function of the 
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three decay angles. As discussed in Sec. 3, the symmetries of the distribution can then be 
used to reduce the number of unknowns; if we consider the real and imaginary amplitude 
components separately, the twelve parameters can be reduced to eight using the symmetry 
constraints. A further spin-amplitude component may be removed by noting that Eq. (2.1) 
is only sensitive to relative normalisations. This leaves seven free parameters at each point 
in q^. However, in [6], only three, out of four, symmetry constraints were considered 
meaning that, in principle, the fits presented were under-constrained. The implications of 
this will be investigated in this section. 

Despite the large increases in B^ — )• K*^fi~^fi~ statistics expected at LHC6, the number 
of signal events available will still be too small for a fixed q^ approach to be taken. Instead, 
the spin-amplitude components are parametrised as second-order polynomials in the region 
q^ E [1, 6] GeV^. These are normalised relative to the value of Re(^o ) ^^ ^ fixed value, Xq, 
of q^. Rather than fitting directly for the amplitudes, we aim to extract the coefficients 
of these polynomials. This introduces a number of model biases: the underlying spin 
amplitudes are assumed to the smoothly varying in the q^ window considered. As noted 
in [6], this was verified for a number of NP models. There is also an implicit assumption 
that the g^-dependent shape of the spin amplitudes is invariant under the symmetries 
of the angular distribution. Neglecting background parameters, the g^-dependent fit has 
((12 — 4) X 3) — 1 = 23 free parameters to be extracted, or 26 in [6]. These will be labelled 
the four- and three-symmetry fits respectively. 

The three-symmetry fit although, in principle, under-constrained is able to converge 
due to the polynomial parametrisation employed. By requiring that three of the spin 
amplitude components vanish for all values of g^, we have used our freedom to choose 
values of 0^, 4>r, and 9 from Eq. (3.10) at each point in g^; the value of 9 is still free to 
vary. However, the PDF, Eq. (2.1), is invariant under changes of 9; hence, the negative 
log-likelihood (NLL) used during minimisation should not be sensitive to its value. The q^ 
dependent shape of each amplitude component is manifestly not invariant under changes 
m 9 - the rotation it implies mixes the imaginary parts of the left- and right-handed 
amplitudes. The polynomial parametrisation of the spin-amplitude components requires 
that each amplitude must be smoothly varying. The fit then selects the value of 9 for each 
signal event which produces the most polynomial-like distribution, as this will have the 
smallest NLL. The general minimising algorithm employed is then able to find a genuine 
minimum and converge properly; the imposition of the polynomial ansatz allowed the 
under-constrained fit of [6] to converge properly. As the experimental observables are 
invariant under all four symmetries, their q^ dependent distributions can be found correctly; 
there are no significant biases seen in the central values extracted compared to the input 
distribution. Small biases are seen in the individual spin-amplitude components; with 
hindsight, correlations between these components were induced by the presence of the 
fourth symmetry. 

4.1 Experimental Analysis 

The discussion above explains why the three-symmetry fit is able to converge successfully, 
and suggests that there should be no change in the experimental uncertainties found when 
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the extra symmetry constraint is added. It is important to demonstrate that this is the 
case. As before, the experimental sensitivity to different observables can be estimated using 
a toy Monte Carlo (MC) approach and used to compare the three- and four-symmetry fits. 

4.1.1 Generation 

An ensemble of data sets for B^ — )■ K*'^n^fi~ can be generated; each data set contains 
the Poisson fluctuated number of signal and background events expected after LHC6 has 
collected 10 fb~ of integrated luminosity. Estimates of the signal and background yields 
were taken from [17, 18] and scaled linearly. The signal distribution was generated using the 
K*^ spin amplitudes discussed in Sec. 2 as input. The contribution from terms including 
the muon mass were included. No assumption of polynomial variation of the amplitudes 
was used in the generation. The signal is assumed to have a Gaussian distribution in nis 
with a width of 14 MeV in a window of niB ± 50 MeV and a Breit-Wigner in mKir with 
width 48 MeV in a window of m^^.o it 100 MeV. A simplified background model is included; 
it is flat in all decay angles, effectively treating all background as combinatorial, but follows 
the q'^ distribution of the signal. Detector acceptance effects as described in [17] are not 
taken into account. When considering CP-conserving quantities, the B and B samples 
are simply considered together. We do not include any contributions from non-resonant 
Bd -^ K~Tr'^fi~^fi~. 

4.1.2 Observable sensitivities 

The ensemble of simulated data sets can then be used to estimate the experimental un- 
certainties expected for a given integrated luminosity at LHC6. For each data set, the full 
angular fit was performed to find the most likely value for each of the free parameters for 
that data set. For the three-symmetry fit there were 27 free parameters; 26 for the signal 
distribution and one to describe the level of background seen. For the four-symmetry fit, 
only 24 parameters were required. In total we created an ensemble of 1200 experiments 
and will, thus, at a given value of g^, get 1200 different determinations of each observable. 
By looking at the point where 33% and 47.5% of results lie within either side of the median 
of the results we can form asymmetric la and 2a errors. Connecting these at different q^ 
values gives us la and 2a bands for the experimental errors on the observable. 

4.1.3 CP asymmetries 

The sensitivity to various CP asymmetries was also considered. In this case, separate B 
and B samples were generated and fit independently. Each sample had on average half the 
number of signal and background events as those described in Sec. 4.1.1. The results of a 
B and a -B fit could then be combined by re-normalising the B amplitudes found, so that 
the extracted value of Re(^Q ) at Xq was the same in both samples. This gives sensitivity 
to CP asymmetries relative to this point. By considering many B and B samples together, 
estimates of the experimental sensitivity to the CP asymmetries could then be found. 
In a real measurement, a more sophisticated approach would be taken which considered 
the two samples simultaneously; however, our simplified approach gives a reasonable first 
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estimate of the experimental sensitivities obtainable and allow comparison with theoretical 
requirements. 

4.2 The polynomial ansatz re-examined 

A key assumption of the fitting approach taken in [6] is that the spin-amplitude compo- 
nents are smoothly varying functions in the range q^ € [l,6]GeV^. It was found that 
when all four symmetries of the massless angular distribution are taken into account, this 
assumption no longer holds; indeed the shape of the spin-amplitude components is not 
invariant under the four symmetries and their shape can be distorted so they are no longer 
well described by second-order polynomials. Other parameterization choices are likely to 
be equally vulnerable to these problems unless they are explicitly invariant under all sym- 
metries of the angular distribution. Consider the three-symmetry case at a fixed q^ value: 
in [6], Aq is removed by setting 9 = arctan(— Ag /j4q) once their phases have been rotated 
away. This can be understood by substituting the trigonometric identities, 

1 

sin(arctan(6l)) = , cos(arctan(6l)) = , (4.1) 

into Eq. (3.10). This introduces a [1 + {Aq /Aq )]~2 term into each non-zero amplitude 
component, which will not be well behaved as Aq — t- 0. For the three-symmetry fit, these 
problems can be avoided by taking Re(^Q ) as the reference amplitude component, forcing 
it to be relatively large at Xq. However, to include the fourth symmetry constraint, a more 
complicated form must be used in order to set four amplitude components simultaneously. 
A different value of each of the four rotation angles is required for every point in g^ due 
to the changing spin amplitudes. There is no guarantee that a set of rotation angles can 
be found such that the unfixed spin-amplitude components resemble smoothly varying 
polynomials for all g^. The q'^ dependence of the SM input amplitude Re(AQ) is shown 
in Fig. 1 once the four symmetries have been applied to fix Im(y4j| ), Im(y4j| ), Re(^j| ), 
and Im(A^) to zero, as required for in the next section. This particular feature is caused 
by Re(An') — ?• at g^ ~ 2GeV^; other rotation choices lead to similar features. The 
distribution can no longer be well described by a second-order polynomial. It may be 
possible to find a choice of rotation parameters that preserve the polynomial features of 
the input spin-amplitude components, however, there are no guarantee that a particular 
choice would work when faced with experimental data. Indeed, an incorrect choice will 
lead to biases in the case where the parametrisation is a poor match for the underlying 
amplitudes. A more generic solution is required and could form the basis for further 
investigations. 

4.3 Fit Quality 

The effect of adding the fourth symmetry constraint was tested, by comparing ensembles 
of three- and four-symmetry fits. The two ensembles were generated with the same random 
seed values so that the ensemble of input data sets was the same for the two approaches. 
The fixed spin-amplitude components were chosen to be lm.{Au), Im(Aj| ), Ile{Ai^), and in 
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the case of the four-symmetry fit also Im(^^). The ampHtudes were stiU normahsed relative 
to Re(^n') at Xq = 3.5 GeV^, however the fits were performed in the range q^ G [2.5, 6] GeV^ 
to avoid the non-polynomial features seen in the spin-amplitude components, such as shown 
in Fig. 1. 

The sensitivities found for the angular 
observables are poorer than those presented 
in [6], due to the decreased signal statistics 
in the reduced q^ window, however, it is 
interesting to compare the performance of 
the two fitting methods. A histogram of the 
NLL of each fit is shown in Fig. 2. The en- 
semble of three-symmetry fits (hatched) and 
four-symmetry fits (solid) can be seen. The 
ensemble of input data sets is slightly dif- 
ferent in each case due to a small number of 
failed computing jobs, but the output distri- 
butions look very similar. This shows that 
the depth of the minima found is approx- 
imately the same for the three- and four- 
symmetry fits. We can also introduce a global correlation factor Gc, which is the unsigned 
mean of the individual global correlation coefficients calculated from the full covariance ma- 
trix. It takes values in the range Gc G [0, 1], where zero shows all variables as completely 
uncorrelated, and one shows total fit correlation. It can be seen in Fig. 3 that the mean 
correlation of the fit is reduced once the fourth symmetry is taken into account. There 
are less outliers at very low Gc and the distribution appears more Gaussian, indicating 
an increase in fit stability has been achieved. The convergence of the fit starting from 
arbitrary initial parameters has also much improved. 



q= (GeV') 



Figure 1: The q^ dependence of Im(AQ ) after 
using the four symmetries of the full-angular dis- 
tribution to fix Im(Ajj'), Im(Aj?), Re(A|), and 
Iin{A^) to zero. 
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Figure 2: The negative log-likelihood factor 
for the three-symmetry (blue hatched) and 
four-symmetry (red solid) ensembles of fits to 
10 fb^^ toy data sets of LHC6 data, assuming 
the SM and with q^ € [2.5,6] GeV^. 



Figure 3: The global correlation factor for 
the three-symmetry and four-symmetry en- 
sembles of fits to 10 fb^ toy data sets of 
LHC6 data, assuming the SM and with q'^ G 
[2.5,6] GeV^. The colour scheme is the same 
as in Fig. 2. 
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Fig. 4 shows the estimated experimental sensitivities found for the theoretically clean 
observable A!^ in the range g^ G [2.5,6]GeV^, with and without the fourth symmetry 
constraint. The fits are for 10fb~ of LHC6 integrated luminosity assuming the SM. As 
might be expected from Fig. 2, there is little difference in the estimated experimental 
resolutions seen. The same conclusion is reached when inspecting other observables. 
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Figure 4: One and two a contours of estimated experimental sensitivity to the theoretically clean 
observable Alj. with full-angular fit to 10 fb of LHC6 data assuming the SM. The results of the 
three-symmetry fit are shown on the left, and the four-symmetry fit on the right. The fits were 
performed in the range q^ G [2.5, 6] GeV^. 



4.4 Discussion 

The discovery of a fourth symmetry in the massless leptons limit of the full-angular dis- 
tribution of B(i —7- K*^ fi^ ii~ requires that the experimental analysis proposed in [6] be 
re-evaluated. The previous analysis used three of the four available symmetry constraints 
to perform a fit, that was, in principle, under-constrained, by parametrising the real and 
imaginary parts of the K*^ spin amplitudes as second-order polynomials. The invariance 
of the observables under all four symmetries, and the freedom to take arbitrary values of 
the 9 rotation angle, allowed the fits to converge and produce correct output, but intro- 
duced a subtle parametrisation bias. As the observables are hy definition invariant to all 
the symmetries, the estimated experimental sensitivities are the same for the two methods. 
This has been demonstrated in this section. However, the need for the development of a 
new fitting method, so that the full experimental statistics available in q^ G [1, 6] GeV^ can 
be used, is now clear. The sensitivities found will be similar to those estimated in Ref. [6] 
and in this paper for the CP asymmetries, but with improved fit stability. 



5. Analysis of CP- violating observables 

In [10, 19], it was shown that eight CP- violating observables can be constructed by com- 
bining the differential decay rates of dT{Bd -> F*°^+£-) and dT{Bd -^ K*° £+£-). In this 
section we analyse the theoretical and experimental uncertainties of those observables in 
order to judge the NP sensitivity of such CP-violating observables. 
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5.1 Preliminaries 

The corresponding decay rate for the CP-conjugated decay mode Bd — )• K*^i~^i~ is given 

by 



d^r 9 



ag^ a cos t// a cos Ok dcp Szvr 

As shown in [10], the corresponding functions Ji{q'^,9i,6K,(p) are connected to functions 
Ji in the foUowing way: 

•^1,2,3,4,7 — ^ «^1,2, 3,4,7, </5,6,8,9 " ^ " «^5,6,8,9 , (5.2) 

where Jj equals J, with ah weak phases conjugated. 

Besides the CP asymmetry in the dilepton mass distribution, there are several CP- 
violating observables in the angular distribution. The latter are sensitive to CP-violating 
effects as differences between the angular coefficient functions, Ji — Ji. As was discussed 
in [10, 19], and more recently in [9], those CP asymmetries are all very small in the SM; 
they originate from the small CP-violating imaginary part of A^ = (VubV*^) / {VtbV^*^) . This 
weak phase present in the Wilson coefficient Cg is doubly-Cabibbo suppressed and further 
suppressed by the ratio of the Wilson coefficients (3Ci +C2)/Cq « 0.085. 

Moreover, it is important to note [19, 9] that the CP asymmetries corresponding to 
J7_8,9 are odd under the transformation (p — )• —cp and thus, these asymmetries are T-odd (T 
transformation reverses all particle momenta and particle spins) while the other angular 
CP asymmetries are T-even. T-odd CP asymmetries are favoured because they involve 
the combination cos{69) sm{6(pw) of the strong and weak phase differences [19, 9], thus, 
they are still large in spite of small strong phases as predicted for example within the 
QCDF/SCET approach. In contrast, T-even CP asymmetries involve the combination 
sin((56') cos(<50vy) [19, 9]". 

Another remark is that the CP asymmetries related to ^5,6,8,9 can be extracted from 
(dr -|- dT) due to the property Eq. (5.2), and thus can be determined for an untagged 
equal mixture of B and B^, mesons. This is important for the decay modes B'^ — t- K*^(-^ 
K^7r^)£~^£~ and Bg — )• <p{-^ K~^K~)£~^£~ but it is less relevant for the self-tagging mode 
Bd -^ K*°{-^ K+ IT- )£+£-. 

Recently, a QCDf/SCET analysis of the angular CP-violating observables, based on 
the NLO results in [13, 14], was presented for the first time [9]. The NLO corrections are 
shown to be sizable. The crucial impact of the NLO analysis is that the scale dependence 
gets reduced to the 10% level for most of the CP asymmetries. However, for some of them, 
which essentially start with a nontrivial NLO contribution, there is a significantly larger 
scale dependence. The g^-integrated SM predictions are all shown to be below the 10~^ 
level due to the small weak phase as mentioned above. The uncertainties due to the form 
factors, the scale dependence, and the uncertainty due to CKM parameters are identified 
as the main sources of SM errors [9] . 



" We note here that this specific behaviour of T-odd and T-even observables was shown in many examples 
of T-odd CP asymmetries (see [20] and references therein) but a general proof of this statement is still 
missing to our knowledge. 
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5.2 Phenomenological analysis 

The NP sensitivity of CP-violating observables in the mode Bd — )• K*^£^i~ was discussed 
in a model-independent way [9] and also in various popular concrete NP models [7] . It was 
found that the NP contributions to the phases of the Wilson coefficients Cy, Cg, and Cio 
and of their chiral counterparts drastically enhance such CP-violating observables, while 
presently most of those phases are very weakly constrained. It was claimed that these 
observables offer clean signals of NP contributions. 

However, the NP reach of such observables can only be judged with a complete analysis 
of the theoretical and experimental uncertainties. To the very detailed analyses in [9, 7] 
we add the following points: 

• We redefine the various CP asymmetries following the general method presented in 
our previous paper [6] : an appropriate normalisation of the CP asymmetries almost 
eliminates any uncertainties due to the soft form factors which is one of the major 
sources of errors in the SM prediction. 

• We explore the effect of the possible AQCB/Tnt corrections and make the uncertainty 
due to those unknown AQCD/^TT-b corrections manifest in our analysis within the SM 
and NP scenarios. 

• We investigate the experimental sensitivity of the angular CP asymmetries using a 
toy Monte Carlo model and estimate the statistical uncertainty of the observables 
with statistics corresponding to five years of nominal running at LHCb (10 fb" ) using 
a full angular fit method. 

We discuss these issues by example of the two angular asymmetries corresponding to the 
angular coefficient functions Jqs and Jg; 

'' ~ d{T + T)/dq^ ' ' ~ d{T + T)/dq^ ' ^^"^^ 

Within the SM the first CP asymmetry related to Jgs turns out to be the well-known 
forward-backward CP asymmetry which was proposed in [21, 22]. 

As a first step we redefine the two CP observables. We make sure that the form factor 
dependence cancels out at the LO level by using an appropriate normalisation: 

j^V2s ^ Jes - J_6s ^ ^v^Js^zjIi, (5.4) 

J2s + J2s Js + J^ 

The Ji are bilinear in the K* spin amplitudes, so it is clear from the LO formulae Eq. (2.6) 
that, following the strategy of [6], any form factor dependence at this order cancels out 
in both observables. We note that J2s has the same form factor dependence as Jgs but 
has larger absolute values over the dilepton mass spectrum that stabilises the quantity. In 
Fig. 5 the uncertainty due to the form factor dependence is estimated in a conservative way 
(see Appendix B) for Aqs defined in Eq. (5.3) and for A^^ defined in Eq. (5.4). Comparing 
the plots, one sees that with the appropriate normalisation, this main source of hadronic 
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Figure 5: SM prediction of the CP-violating observables Aqs (left) and A^^'^ (right) as function 
of the squared lepton mass with uncertainty due to the soft form factors only. Notice the difference 
in scale and the difference in relative error in the two figures. 




3 4 

Figure 6: SM prediction of the CP-violating observables A^ (left) and Ag (right) with uncertainty 
due to the soft form factors only. Notice the difference in scale of the two figures. 

uncertainties gets almost eliminated. The leftover uncertainty enters through the form 
factor dependence of the NLO contribution. Fig. 6 shows the analogous results for the 
observable ^g . 

In the second step we make the possible Aqcd/"^6 corrections manifest in our final 
results by using the procedure described in Sec. 2.3. It turns out that in spite of this very 
conservative ansatz for the possible power corrections (we neglect for example any kind of 
correlations between such corrections in the various spin amplitudes), the impact of those 
corrections is smaller than the SM uncertainty in case of the two observables Aq^ and Ag . 
In the left plot of Fig. 7 the SM error is given, including uncertainties due to the scale 
dependence and input parameters and the spurious error due to the form factors. In the 
right plot the estimated power corrections are given, which in case of the CP-violating 
observable Aq^ are significantly smaller than the combined uncertainty due to scale and 
input parameters. Fig. 8 shows the same feature for the CP-violating observable Ag . This 
result is in contrast to the one for CP-averaged angular observables discussed in [6], where 
the estimated power corrections always represent the dominant error. The reason for this 
specific feature is the smallness of the weak phase in the SM. Thus, one expects that the 
impact of power corrections will be significantly larger when NP scenarios with new CP 
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Figure 7: SM uncertainty in Aq^'' (left) and estimate of uncertainty due to vlgcD/'Tih corrections 
with Ci^2 = 10% (right). 
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Figure 8: SM uncertainty in A^ (left) and estimate of uncertainty due to vIqcd/w;, corrections 
(right, light grey (green) corresponds to Ci,2 = 5%, dark grey (green) to Ci,2 = 10%). 

phases are considered (see below). 

In the third step we consider various NP scenarios. Here we follow the model-independent 
constraints derived in [9] assuming only one NP Wilson coefficient being nonzero. We con- 
sider three different NP benchmarks scenarios of this kind: 



1. 



|C^P| = 2 and (/.^P 



!' 2 



, TT (Red) ; 



n |/^NP| 



1.5 and (/)5^o'' = f>f.^ (Grey); 
3. |C;^g| = 3 and (p^Q = f ) f )^ (Blue); 

where the colours refer to the ones used in the following figures. The absolute values 
of the Wilson coefficients are chosen in such a way that the model-independent analysis, 
assuming one nontrivial NP Wilson coefficient acting at a time, does not give any bound 
on the corresponding NP phase. 

Fig. 9 shows our two observables in the three scenarios with the phase value |: the 
CP- violating observable Aq^ might separate a NP scenario (2), while the central values of 
scenarios (1) and (3) are very close to the SM. Moreover observable ^g' seems to be suited 
to separate scenarios (1) and (3) from the SM. 

However, to judge the NP reach we need a complete error analysis within the three 
NP scenarios. As shown in Sec. 2.3 we now work with three weak sub-amplitudes in which 
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Figure 9: NP scenarios, assuming one nontrivial NP Wilson coefficient at a time, next to SM 
prediction for A^^'^ (left) and A^ (right), for concrete values see text. 

possible power corrections are varied independently. The plots in Figs. 10 and 11 show that 
the possible ^QCD/^^-fe corrections have a much larger impact on our two observables in 
the NP scenarios than in the SM and become the dominating theoretical uncertainty. We 
also get significantly larger possible ^Qco/w-b corrections when changing the value of the 
new weak phase from f to |. Regarding even larger phase values, we note here that the 
NP effects drastically decrease again when phase values around tt are chosen as expected. 
Nevertheless, in view of the theoretical ^QCD/'Tift uncertainties only, the two CP-violating 
observables could discriminate some specific NP scenarios with new CP phase of order | 
or ^ from the SM; in case of ^g/* NP scenario 2, in case of Ag NP scenario 3 and possibly 
1. 




q'(GeV') 



q'(Ge\f) 



Figure 10: A^^'^: Estimate of uncertainty due to AQcu/mi, corrections within NP scenarios as in 
previous figure with phases 0i = f (left) and (f>i = ^ (right). 

One should also consider the additional theoretical uncertainties due to scale depen- 
dence, input parameters and soft form factor dependencies within the NP scenarios. Those 
additional theoretical uncertainties are sizable and of the same order as the ones due to 
^QCD/'iT'b corrections: they are shown in the left plots in Figs. 12 and 13 as orange bands 
overlaying the total errors bars including also the ylQco/'Tib corrections. 

As the last step, we analyse the experimental sensitivity of the angular CP asymmetries 
using a toy Monte Carlo model. The right plots in Figs. 12 and 13 show the estimates of 
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Figure 11: A^: Estimate of uncertainty due to AqcB/'mb corrections within NP scenarios as in 
previous figure with phases 0i = f (left) and (t>i — ^ (right). 




12 3 4 5 6 

cfiGeV') 

Figure 12: Aq^^: Estimate of uncertainty due to ^QCD/^JT-b corrections (grey bands) in NP scenario 
2, \Ciq\ = 1.5 and (p^Q = | with the other theoretical uncertainties overlaid (orange bands) and in 
SM (left) and experimental uncertainty (right). 



the statistical uncertainty of Aq^ and Ag with statistics corresponding to five years of 
nominal running at LHCb (10fb~ ). The inner and outer bands correspond to la and 2a 
statistical errors. The plots show that all the NP benchmarks are within the la range of the 
expected experimental error in case of the observable Aq^, and within the 2a range of the 
experimental error in case of the observable A^ . We emphasise that from the experimental 
point of view the normalisation is not important when calculating the overall significance 
because the overall error is dominated by the error on the numerator. So the experimental 
error of the observables Aqs and A^ defined in Eq. (5.3) using the traditional normalisation 
will be similarly large to the one of our new observables Aq^ and Ag defined in Eq. (5.4). 

Our final conclusion is that the possibility to disentangle different NP scenarios for 
the CP- violating observables remains rather difficult. For the rare decay Bd — t- K*^£~^£~, 
LHC6 has no real sensitivity for NP phases up to values of | (and neither up to values 
of vr) in the Wilson coefficients Cg, Ciq and their chiral counterparts. Even Super-LHC6 
with 100 fb~ integrated luminosity does not improve the situation significantly. This is 
in contrast to the CP-conserving observables presented in [6] and further discussed in the 
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Figure 13: A^ : Estimate of uncertainty due to ^qcd/"^6 corrections in NP scenarios 1 ( \C, 



ri 



br 



§, red bands) and 3 (|C 



j^gl =3, (j)iQ = ^, blue bands) with the other theoretical uncertainties 
overlaid (orange bands) and in SM (left) and experimental uncertainty (right). 

next chapter which, both from the theoretical and experimental point of view, are very 
promising. 

6. Analysis of CP-conserving observables 

The CP-conserving observables can be analysed at LO in the large recoil limit using the 
heavy-quark and large- i?;^'* expressions for the spin amplitudes, as first proposed in [4]. 
One of the advantages of this approach is that we obtain analytic expressions of these 
observables in a very simple way. These expressions can be used to study the behaviour of 
the observables without having to rely on numerical computations, since the most relevant 
features arise already at LO. The main goal of this section is to perform this type of analysis 

(i) 

on the ^rp'^ observables. 

(2) 

6.1 Leading-order expressions of A^ 

(2) 

The asymmetry Al^ , first proposed in [4] is given by 

l^iP-l^iiP 



A, 



(2) 



M,|2 + |A||I2^ 



H 



(6.1) 



where \Ai\'^ = \A^\'^ + \A^\'^. It has a simple form, free from ^_l(0) form factor dependencies, 
in the heavy-quark (m^ — )■ oo) and large K*'^ energy {Ex* — )• oo) limits**: 



,(2) 



Re CioCi*o + P'Re C^q + FRe CjC^ 



^ |Cio|2 + ICioP + P2(|C7|2 + |C;|2) + ICsP + 2PRe (C7C9*) ' 
where F = 2rrnymB/q^- The Wilson coefficients can take the most general form: 



Ci = Cf" + \C, 



NP|gt(^ 



C' = \C'W 



7,9,10. 



(6.2) 



(6.3) 



"Notice that along this section we will drop the superscript "efF' that C7 and Cg should bear in order to 
simplify the notation. 
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We will neglect henceforward both the tiny SM weak phase (pg^, that arises from the 
CKM elements ratio Xu = (KbK's)/(^bKs)) ^^d the SM strong phase 6*9'^, smaller than 
P in the low dilepton mass region 1 GeV^ ^ q'^ ^ 6 GeV^ [22]. 

Obviously, the observable AJp vanishes in the heavy-quark and large K*^ energy lim- 
its at LO when all the Wilson coefficients are taken to be SM-like. This result can be 
understood rather easily. The left-handed structure of weak interactions in the SM guar- 
antees that, in these limits, the s quark created in the b ^ s transition will have helicity 
h{s) = —1/2 in the massless limit {nis — )• 0) [23]. This s quark will combine with the 
spectator quark d of the Bd to form the K*^ meson with h{K*^) = — 1 or (but not +1), 
therefore Hj^ = at quark level in the SM. Using Eq. (2.3), this translates into A±_ = —Ah 
at the quark level, which corresponds to A±_ ~ ~^\\ at the hadron level [24, 25, 26]. 

The NP dependence of Ar^ can be stud- 
ied in a model independent way by switch- 
ing on one Wilson coefficient each time and 
keeping all the others at their SM values. 
A simple inspection of Eq. (6.2) shows that 
only the chirally flipped operators Oj and 
OiQ give a non-zero expression for AJj< in 
our approximation: 

2F{FCf^ +Cf')\C'rj\cos{4>'j) 



1.0 



0.5 



0.0 



-0.5 



-1.0 




a(2) 

and 

.(2) 



(Cf^)^ + F^lCjl"^ + {FC 



q'{Ge'/) 



SM _j_ /'SM\2 

(6.4) 



.(2) 



n /-"SMI/i' 

^•-10 I'-io 



COS 



^lOJ 



10' 



(Cf5*)2 + |C;o|2 + (FCf^+C9S^)2 

(6.5) 



(2) 



Equations. (6.5) and (6.4) show that A)^ is 
sensitive to both the modulus and the sign 
of the Wilson coefficients Cj and Ciq. When 
NP enters only C^q, the fact that Cio < 
in the SM makes the observable negative 
< IT, enabling us to distin- 



unless I < 



^101 



Figure 14: A^ ' in the SM (green) and with NP 
in C^Q = 3e't (blue), this value is allowed by the 
model independent analysis of [9]. The inner line 
corresponds to the central value of each curve. 
The dark orange bands surrounding it are the 
NLO results including all uncertainties (except for 
^QCD/"^b) as explained in the text. Internal light 
green/blue bands (barely visible) include the esti- 
mated v1qcd/"t-& uncertainty at a ±5% level and 
the external dark green/blue bands correspond to 
a ±10% correction for each spin amplitude. 



,(2) 



guish the sign of this weak phase (Fig. 14). Likewise, if NP appears in C7, A)^' will display 
a zero in the dilepton mass spectrum when FC^^ -\- Cg"^ = 0, which will coincide exactly 
with the zero of the observable A-p-Q at LO [13]. As the zero is independent of Cy, all curves 
with Cf^ should exhibit it at q^ ~ 4GeV2, but if there is also a NP contribution to C7, the 
zero will be shifted either to higher or lower values of cf'. In case of a sign flip affecting 
C7, A)^ would not have a zero at any value of q^, exactly as for ApB (see [27] for a recent 
discussion of different mechanisms to achieve this). In fact, should NP enter both Oj and 
O-j simultaneously, Eq. (6.2) would imply 



A. 



(2) 



J I 7NP 



oc2F 



(Fcr+cf") 



C7|cos( 



+ F\C 



|Cricos( 



') 



(6.6) 



28 



while 



A 



FB 



■jl yNP 



o^ FCf + Cf^ + F\C\ 



cos 



^T) 



(6.7) 



The comparison of Eq. (6.6) with Eq. (6.7) can be used to explain the improved sensitivity 



i(2) 



of A)j. to certain types of NP versus that of Afb- The numerator of A)j^ exhibits sensitivity 



(2) 



to the weak phases (p^^ and 



'7^ 



having an interference term enhanced by the large factor F 



(8 ^ -^ ^ 48 in the dilepton mass region studied), while ^fb is only sensitive to i;^^''. Thus, 
a wider departure from the SM behaviour is to be expected in A^ when NP enters the 
operators 0^ and O-j. This is shown in Fig. 15 using three different scenarios, described in 
the caption of Fig. 15, compatible with present experimental and theoretical constraints. 
Therefore, we emphasise that A)j, must be regarded as an improved version of ^fb once 
the full-angular analysis becomes possible. 
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3 4 



Figure 15: Observables ^r^ and Afb with NP curves for three allowed combinations of C7 and 
C7 following the model independent analysis of [9]. The bands correspond to the SM and the 
theoretical uncertainty as described in Fig 14. The cyan line (shown with the label a) corresponds 
to (C^^, C7) = (0.26e-'TS, 0.2e*^), the brown line h to (0.07e*T, 0.3e'T) and the magenta line c 
to (O.OSe*'^, 0.07). 



(5) 

6.2 Leading-order expressions oi A>r^ 

In the SM, we get in the heavy-quark and large- -Ex* limits at LO: 



A 



(5) 



-(Cio 



SM\2 



+ (FCf * 



_1_ /'SM\2| 



SM 



2[(Cf^)2 + (FCf^+qM)2]' 



,(5) 



which sets the "wave-like" behaviour of Ai^ . At low g , Eq. (6.8) can be used to check 



that A 



(5) 1 1 GeV^ 



(2) 



,(5) 



0.4. On the other hand, at the zero-point of A)^ and Afbj ^t exhibits 

,(5)|4GeV2 



an absolute maximum of magnitude A>r^ V^." ~ 0.5. 

Any inclusion of NP in the Wilson coefficients C7, Cg and Cio will give rise to the 
appearance of an extra term in the numerator (with respect to Eq. (6.8)) that will shift 
the observable along the y-axis. 



A 



(5) 



7r/2 
7NP 



-(cf^Y + F2|CNP|2 + {pcf' + c^^Yf + 4 {FC^\c. 



NP|]2 



2 [(Cf5')2 + F2|C^P|2 + [FCf^ + Cl^Y\ 



(6.9a) 
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Figure 16: A)^' in the SM and with NP in C^q = 3e's and C^'' = 2e'8 (left) and in both C7 and 
C7 Wilson coefficients (right). The cyan hne (a) corresponds to (Cy^, C-j) = (0.26e~'TO, 0.2e"), the 
brown hne (6) to (0.07e*T, 0.3e*T) and the pink hne [d) to (0.18e~*5, 0). The bands symbohse 
the theoretical uncertainty as described in Fig. 14. 



A, 



(5) 



A 



(5) 



7r/2 
9NP 

7r/2 
IQNP 



-(Q02 + |CNP|2 + (^CfM + qM)2]2 + 4 [CSM|CgNP|]2 



/hc 



2[(Cf^)2 + |qp|2 + (FCf 



.CSM)2] 



SMN2 

10 / 



ICfol^ + {FCf^ + C9s^)2]2 + 4 [|Cfon(^Cf" + C, 



SM)]2 



(6.9b) 



(6.9c) 



2 [(Ci^^)2 + ICfJ-p + (FCf^ + C9^^^)2] 

In Eq. (6.9) we have chosen for simplicity the weak phase (f)^^ = 7r/2 for z = 7, 9, 10, but 
they turn out to be dominated by the SM contribution unless the NP Wilson coefficients are 
very large. However, if the weak phases associated to NP Wilson coefficients are different 
from 7r/2, the A)j. curve will get shifted either to the left or to the right, depending on the 
value of the angle, as shown in Fig. 16. 

NP might also enter via the chirally flipped Oj and O^g. The corresponding LO 

(5) 

expressions of A)^, in the heavy-quark and high-i?/<* limits read 

-{Cl^f + iFq^ + Cf'f 



A. 



(5) 



f2|c;|2 



2 [{ci^Y + {Fq^^ + CSM)2 + F2|c;|2 



and 



A. 



(5) 



-(Cl^y + \C,,\^ + {FCr+Cr) 



SM\2 



10' 



2[icf^r + \c[,\^ + {Fcr+ci 



SM^2 



(6.10) 



5.11) 



(5) 

Equations (6.10) and (6.11) are both free from NP weak-phase dependence, ^rp"^ evaluated 
at the g2 value of the ^pB zero-point can be computed easily using Eq. (6.11), obtaining 



A. 



(5)1 



"?g 



/'/iSM\2 I \p' |2| 
~ I'-IO ) + I'-'IOI I 

(c-)2 + |c;oP 



(6.12) 



where the choice C^q = enables us to recover the SM prediction ^L L,!* = 0.5. In 



Fig. 16 (left) it can be seen that for |C^q| = 3 the departure of the NP curve obtained from 
the SM behaviour is indeed large. 
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6.3 Analysis of a'^^ and A^^^ 

The observables Al^ and Arj,' were first introduced in [6] to test the longitudinal spin 
ampUtude Aq in a controlled way: 

Unfortunately, the simultaneous appearance of A±, A\i and Aq inside square roots turns 
the heavy-quark and large-energy limits into rather awkward expressions, not really useful 
to explain the behaviour of these observables at a glance. Therefore, we only outline 
their general properties. Equation (6.13) shows that A!^ and A)^' play a complementary 
role, as the numerator of A^ and the denominator of A^ are the same. Thus, when a 
minimum appears in one of them, a maximum is expected in the other observable and the 
other way around. This is indeed what can be observed in Fig. 17. For the values of the 
Wilson coefficients chosen, NP entering C^g can easily be distinguished from the SM curve, 
displaying a maximum at around 3.5-4 GeV^ (exactly in the energy region where A^ 
is showing a minimum), while Ciq can only be clearly identified using Al^ . Something 
similar happens with NP entering Cj^ and Cj: the model-independent values chosen for 

(3) 

these Wilson coefficients do not give rise to clear NP signals from A}^, , but they can be 
easily told apart using A!^. . In those situations where the origin of the NP curve can not 
be clearly established using a single observable (for instance, the c curve in the A^ plot 
of Fig. 17 is very similar to the C^q curve), the combined use of A^ , Al^ , A)^ , A^p and 
maybe ^fb enables us to identify which Wilson coefficient (s) has a contribution from NP. 



7. Conclusion 

In this paper we have presented how the decay B^ — )■ K*^i~^£~ can provide detailed knowl- 
edge of NP effects in the flavour sector. We developed a method for constructing observables 
with specific sensitivity to some types of NP while, at the same time, keeping theoretical 
errors from form factors under control. A method based on infinitesimal symmetries was 
presented which allows in a generic way to identify if an arbitrary combination of spin 
amplitudes is an observable of the angular distribution. For the case of massless leptons 
we identified the explicit form of all four symmetries present. We showed the possible 
impact of the unknown AQCj^/nib corrections on the NP sensitivity of the various angular 
observables in a systematic way using an ensemble method. Experimental sensitivity to the 
observables was evaluated for datasets corresponding to 10fb~ of data at LIIC6. Using 
these tools, we did a phenomenological analysis for both CP-conserving and CP-violating 
observables. The conclusion from this is that the CP-violating observables have very poor 
experimental sensitivity while the CP-conserving observables A)^ (with i = 2, 3, 4) are very 
powerful for finding NP, including situations with large weak phases. 
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Figure 17: A^r^^ and Ai^'> in the SM and with NP in and Cf(f = 1.5e*t and C^q = 3e'i (left) and 
in both C7 and C7 Wilson coefScients (right). The cyan line (curve a) corresponds to {C^^ , C7) = 
(0.26e~'Tlj 0.2e"), the brown line (curve b) to (0.07e*T, O.Sc't) and the magenta line (curve c) 
to (O.OSe*'^, 0.07). The bands symbolise the theoretical uncertainty as described in Fig. 14. 
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A. Kinematics 



Assuming the K*^ to be on the mass sheh, the decay Bd — )■ K*^i'^i~ is completely de- 
scribed by four independent kinematic variables; namely, the square of the lepton-pair 
invariant mass, q'^, and the three angles 6i, 9k and cj). The sign of the angles for the B^ 
decay shows great variation in the literature. Therefore we present here an explicit defi- 
nition of our conventions and point out where the same or different definitions have been 
used. 

First we consider the B^ — > K*^i'^i^ decay. The angle 9i is the angle between the /i"*" 
momentum in the rest frame of the dimuon and the direction of the dimuon in the rest 
frame of the B^. The 6k angle is in a similar way the angle between the K^ momentum 
in the K*^ rest frame and the direction of the K*^ in the rest frame of the Bd- 

Let us for Bd — )■ K*^i'^l~ define the momentum vectors 



Pi+i- = Pe+ +Pe- , (A.l) 

Qi+£- = P£+ - Pi- , (A. 2) 

Pk-7t+ = Pk- +P7T+ , (A.3) 

Qk-7t+ = Pk- -Ptv+ ■ (A.4) 

In the dimuon rest frame, we have that the £"*" momentum is parallel to Qi+^- and also 

that Pk--k+ points in the opposite direction of the dimuon in the Bd rest frame. Thus we 
can compute the 9i angle as 

'~ in« iip« I' ^ ' 

where the superscript is used to indicate the frame. In a similar way we have in the K*^ 
rest frame 

\^K-TT+\\-^i+e.-\ 

Finally, if we go to the rest frame of the Bd, we have (j) as the signed angle between 

the planes defined by the two muons and the K*^ decay products respectively. Vectors 
perpendicular to the decay planes are 

which lets us define (p from 

Ng+i- ■ NK-Tr+ . , I Nf+e- X ^K-TT+\ Pe+e- /a q\ 

cos(j) = = = , sm(j) = \ —^ ^ • _ — . (A. 8) 

\N,+t-\\NK-^+\ \\Nf+,-\\NK-^+\) \P^+,-\ 

The angles are defined in the intervals 

-l^cos6';^l, -li^cosOK <,l, -7r^(/><7r. (A.9) 
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The definition given liere is identical to [6] but is different to [7]. However, the two defini- 
tions result in the same signs for all the coefficients Jj in Eq. (2.4). 

Now for the B^ — t- K*^i~^£~ decay the 9i angle is still specified with respect to the £"*" 
while for 6k the angle is for the K~^. This is equivalent to what is done in [7]. As the 6i 
angle does not change the sign of the lepton, we have 

•^1,2,3,4,7 = Jl,2,3,4:J ' <^5,6,8,9 = " o^5,6,8,9 • (A. 10) 

in the full-angular distribution in the absence of CP violation. 

For the experimental papers [1, 2], a definition has been adopted where all angular 
distributions have been plotted for the B^ — )• K*^£~^i~ decay, with the Bd -)• K*^i~^£~ 
events overlaid assuming CP conservation. In practise this means that B^ — )■ K*^i~^£~ 
events have the sign of cos 6i reversed before plotting. When experiments progress to 
measuring the (p angle as well, special care needs to be taken to get the definitions correct. 

B. Theoretical input parameters and uncertainties 

To compute the soft form factor error bands in Figs. 5 and 6 in a conservative fashion, we 
have used, as input data, the values of 6|(0) and i^_l(0) shown in Table 2. One can notice 
that the ^_l(0) value has been taken from [7], as it is compatible to ■^±(0) = 0.26 used in 
[14], while for ^m (0) we have kept the value from [14] to allow for a wider uncertainty range. 
The g^-dependence of the form factors V, Ai and A2 has been parametrised according 
to [28] 

where i^(0), ap and bp are the fit parameters shown in Table 3 of [28]. Substituting the 
outcomes of Eq. (B.l) into [14] 

rriB + rriK* 

^11 W = — TTB ^19 M{q , (B.2 

2Ek* ms 

we can obtain both the central value and the associated uncertainty curves for 61(9^) and 
C±iQ'^) ill the 1-6 GeV^ range. These are used to get the fitting parameters A, B, C and D 
of 



A-B{qyml) 



2^ ^ ^--^ ' 1 



where A, C ~ 1 within a per mille precision. This parametrisation follows closely Eq. (47) 
in [13] and allows us to explore the impact of ^±(0) and 61 (0) (with their corresponding 
uncertainties) to the CP- violating and CP-conserving observables studied throughout this 
paper. 
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TUB [29] 


5279.50 ± 0.30 MeV 


A [29] 


0.226 ±0.001 


ruK [29] 


896.00 ± 0.25 MeV 


^[29] 


0.814 ±0.022 


Mw [29] 


80.398 ± 0.025 GeV 


p[29] 


0.135 ±0.031 


Mz [29] 


91.1876 ± 0.0021 GeV 


7? [29] 


0.349 ±0.017 


mt{mt) [13] 


167 ± 5 GeV 


4£'^ [29] 


220 ± 40 MeV 


rhb{mb) [30] 


4.20 ±0.04 GeV 


as{Mz) [29] 


0.1176 ±0.0002 


rhci'rhc) [31] 


1.26 ±0.02 GeV 


ftem [29] 


1/137 


/b [32] 


200 ± 25 MeV 


ai'|.(2GeV) [33] 


0.03 ±0.03 


/^.(2GeV) [33] 


163 ± 8 MeV 


a^' II. (2 GeV) [33] 


0.08 ±0.06 


/)l,(2GeV)[33] 


220 ± 5 MeV 






mBe|l(0)/(2mA^.) [14] 


0.47 ±0.09 


AB,+(^ih) [34] 


0.51 ±0.12 GeV 


e±(0) [7] 


0.266 ±0.032 


/^U7] 


2.2 GeV 



Table 2: Summary of input parameters and estimated uncertainties. 

The next step is to compute the amplitudes, keeping one soft form factor fixed at the 
central value and varying the other in the range allowed by its uncertainty. From them, the 
observables can be obtained in a straightforward way and the errors added in quadrature. 

To generate the theoretical error bands not due to A/rrib corrections (plotted as the 
inner orange strips in the plots of Sees. 5 and 6) we have used the criteria of Beneke et 
al. in [13] and added the following uncertainties in quadrature: the renormalisation scale 
uncertainty has been found by varying ^ between 2.3 and 9.2 GeV (where /x is the scale at 
which the Wilson Coefficients, Us and the MS masses are evaluated), the uncertainty in 
the ratio mc/rrib by varying this quantity between 0.29 and 0.31, and the other parametric 
uncertainties have been collected into the factor [6] 



<l' 



NcruBUQ^) 



with z =±, 



(B.4) 



that determines the relative magnitude of the hard-scattering versus the form factor term 
[13], which is uncertain by about ±35%. In our numerical analysis we have used the values 
of the Wilson coefficients in Table 1 of Ref. [13]. 
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